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Abstract 

A previously derived relativistic energy density functional for nuclei, based on 
low-energy in-medium chiral dynamics, is generalized to implement constraints from 
chiral SU (3) effective field theory and applied to A hypernuclei. Density-dependent 
central and spin-orbit mean fields are calculated for a A hyperon using the SU{3) ex- 
tension of in- medium chiral perturbation theory to two- loop order. Long range AA^ 
interactions arise from kaon-exchange and from two-pion-exchange with a S hyperon 
in the intermediate state. Short-distance dynamics is encoded in contact interac- 
tions. They include scalar and vector mean fields reflecting in-medium changes of 
quark condensates, constrained by QCD sum rules. The A single particle orbitals 
are computed for a series of hypernuclei from ^^C to ^'^^Pb. The role of a surface 
(derivative) term is studied. Its strength is found to be compatible with a corre- 
sponding estimate from in-medium chiral perturbation theory. Very good agreement 
with hypernuclear spectroscopic data is achieved. The smallness of the A- nuclear 
spin-orbit interaction finds a natural explanation in terms of an almost complete 
cancellation between short-range scalar / vector contributions and longer range terms 
generated by two-pion exchange. 
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1 Introduction 



The spectroscopy of A hypernuclei [T]|2|5] sets strong constraints [Ij on tlie 
spin dependence of tlie AA^ effective interaction. There is convincing evidence 
that the A-nuclear spin-orbit interaction is abnormally weak compared to the 
very strong spin-orbit force experienced by nucleons in ordinary nuclei. Ex- 
periments setting limits on the spin-orbit coupling in light A hypernuclei are 
reported in Ref. [3]. A recent analysis [B] of heavy hypernuclei confirms the 
systematic smallness of the A-nuclear spin-orbit splitting. 

Many theoretical attempts were made over the years to understand the weak- 
ness of the A-nuclear spin-orbit coupling [TfSlQfTOlITT] . with model assumptions 
ranging from quark degrees of freedom via phenomenological boson exchange 
mechanisms to unusually strong (negative) uoK tensor couplings. The system- 
atics of A-nuclear shell model orbitals have been studied in a variety of ap- 
proaches. The quark-meson coupling model p^][T3] seeks the origin of the weak 
spin-orbit coupling of the A at the quark level. Relativistic mean-field (RMF) 
models [7|[T^ generate the spin-orbit coupling by the coherent interplay of 
scalar and vector mean fields. Even though these fields are only about half 
as strong for A hyperons as compared to those for nucleons in nuclei, the 
resulting A-nuclear spin-orbit splittings are usually still too large unless ad- 
ditional ad-hoc mechanisms are invoked for their suppression. In mean-field 
descriptions with phenomenological density dependent interactions of Skyrme 
type [1^161117] . and in several hypernuclear many-body calculations (Fermi 




Hypernetted Chain, Brueckner-Hartree-Fock) p^SflQ] , the spin-orbit splitting 

is suppressed simply by hand. 

In this article we follow a different path motivated by recent developments at 
the interface between low-energy QCD and the nuclear many-body problem. 
Chiral effective field theory is established as the realization of QCD in the 
low-energy limit, with pions as Goldstone bosons of spontaneously broken 
chiral symmetry being the active light degrees of freedom. In its version with 
nucleons and A(1230) isobars as heavy baryons, chiral effective field theory 
is considered to be an appropriate starting point not only for describing low- 
energy pion-nucleon interactions but also for the nuclear many-body problem 



In-medium chiral perturbation theory [2T|22] provides a successful framework 

for constructing the energy density of nuclear matter as a function of Fermi 
momentum kp. Long and intermediate range one- and two-pion exchange in- 
teractions including tensor and three-body forces are treated explicitly. Short- 
distance dynamics, not resolved in detail at the characteristic Fermi momen- 
tum scales, is encoded in contact terms. The contact terms generate contribu- 
tions to the energy per particle that are linear in the density p = 2kp/ (Svr^). 
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These terms need to be adjusted, e.g. by reproducing the empirical binding 
energy, while terms of higher (fractional) powers in the density, involving ef- 
fects of the filled Fermi sea on two-pion exchange processes, are computed 
with input fixed entirely by empirical low-energy pion-nucleon data. Three- 
nucleon interactions are systematically incorporated. Spin-orbit forces follow 
consistently from the evaluation of spin-dependent terms in inhomogeneous 
nuclear matter [23|24ll25ll26|27j . 

The translation of this framework into a relativistic energy density functional 
for finite nuclei is designed so that it keeps the separation-of-scales 

concepts of the previous nuclear matter calculations. It includes explicit pion 
and two-pion exchange dynamics up to three-loop order in the energy den- 
sity. It also incorporates strong scalar and vector mean fields (equivalent to 
contact interactions and constrained by in- medium QCD sum rules), and an 
additional surface (derivative) term. Applications of this approach to nuclei 
have been quite successful throughout the nuclear chart and naturally suggest 
a systematic extension to hypernuclei, with special focus on the very different 
spin-orbit interactions of nucleons and A hyperons. 

The generalization to chiral SU(3) introduces a well-defined set of couplings 
of the pseudoscalar meson octet to the baryon octet. The longest range AA^ 
interactions are now generated by two-pion exchange with an intermediate 
SA^ and (less importantly) by kaon exchange. The even longer range one-pion 
exchange is excluded in lowest order by the isospin J = of the A. Contact 
interactions are again representing unresolved short- distance dynamics. A first 
study using this in-medium chiral SU(3) dynamics approach was performed for 
^^O as a test case [30]. In the present work these calculations are systematically 
expanded over a large range of hypernuclei from ^^C to ^°^Pb. 

The basic mechanisms that govern spin-orbit interactions in comparison be- 
tween nuclei and hypernuclei are a persistent theme of the present study. Three 
major sources of spin-orbit interactions in nuclear systems can be identified: 

(i) short-distance dynamics of coherently acting scalar and vector mean 
fields; 

(ii) intermediate range spin-orbit forces induced by the pion exchange tensor 
interaction in second order, with Pauli blocking of intermediate nucleon 

states; 

(iii) a three-body spin-orbit interaction of Fujita-Miyazawa type ^31j, pro- 
duced by two-pion exchange with intermediate excitation of a virtual A 
isobar. 

All of these mechanisms are generated within the in-medium chiral dynamics 
framework. The important feature pointed out in Refs. [23|32] and further 
elaborated in Ref. [33], is that mechanism (ii) comes with opposite sign but 
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similar magnitude as compared to mechanisms (i) and (iii). The balance be- 
tween all three mechanisms is shown to account for the large residual spin-orbit 
splitting observed empirically in nuclei. For a A in a hypernucleus, however, the 
three-body Fujita-Miyazawa mechanism (iii) has no analogue simply because 
the A exists only as a single valence particle and there is no hyperon Fermi sea. 
This implies that the short-distance (or scalar-vector) mechanism (i) and the 
intermediate range, second-order tensor force mechanism (ii) largely cancel 
[32] to make a small net spin-orbit splitting for the A-hypernuclear orbitals. 
Testing this scenario over a large set of hypernuclei is one of the primary goals 
of the present work. 



2 Theoretical framework 

2.1 Hypernuclear energy density functional 

A reliable and accurate calculational framework to deal with fermionic many- 
body systems such as nuclei and hypernuclei is the density functional ap- 
proach. Here we start from a relativistic energy density functional for hyper- 
nuclei with a single A hyperon orbiting in the nuclear environment: 

E[P]=E''[P]+EL[P]+EL[P], (1) 
where is the energy of the nuclear core and 

eL = J dM%\M-n-'^ + Mj,]ijj,\%) , (2) 

eL = J rfV{($o|Gf(p) (i^^p) (tAa^a) |<I>o) 

+ ($o|^| 5m(^^)^^(^aV^a)I'^o)} , (3) 

are the additional contributions to the energy involving the hyperon. Here 
|$o) denotes the hypernuclear ground state; iPa{x) and tp{x) are the hyperon 
and nucleon fields, respectively. 

The E[p\ of Eqs. ([I]l3]) represents a generalization to hypernuclei of the nuclear 
energy density functional [28|29] , constrained by basic features of low-energy 
QCD. The nuclear part E'^[p], introduced in Ref. [29], describes the core of 
interacting protons and neutrons in terms of the corresponding isoscalar and 
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isovector densities and current^. The hypernuclear functional, Eq. ([T]), in- 
cludes the A kinetic energy and mass term E^^^ of Eq. ([2]), and the term E(^^ 
of Eq. ([3D which summarizes A-nucleon interactions in the nuclear environ- 
ment. 

Also included in this density functional is a surface term proportional to Dg 
that involves gradients of the isoscalar-scalar nucleon density and the scalar 
hyperon density distributions. Such a term arises naturally in the gradient 
expansion of the energy density functional for a finite system. The leading 
pieces (the first two terms on the r.h.s. of Eq. ([3])) correspond to the lo- 
cal density approximation (LDA). They account for the interaction of a A 
with homogeneous isospin-symmetric nuclear matter taken at the actual local 
density of the nuclear core. The A-hypernucleus is spatially inhomogeneous, 
however, and surface effects are expected to be non-negligible. The gradient 
term (third term in Eq. ([3])) is introduced to account for such corrections 
beyond the LDAE]. The explicit form of this surface term is guided by the 
corresponding part of -E^[p] in the nucleon sector [29] . 



2.2 Effective interaction 

The strength of the effective A-nucleon interaction in the nuclear medium, see 
Eq.([3]), is determined by density-dependent vector and scalar couplings, Gy{p) 
and G^{p), of dimension (length)^. We follow here a strategy analogous to our 
previous calculations for nuclei. Long and intermediate range hypernuclear 
dynamics are governed by chiral two-pion and kaon exchanges in the presence 
of the nuclear Fermi sea. Short distance dynamics is encoded in contact terms 
which generate Hartree type mean-field contributions linear in the density. The 
hyperon self-energies resulting from both long and short range interactions 
are then transcribed into equivalent density dependent couplings, expanded 
in (fractional) powers of the local density p{r). 

Consider first the effective couplings generated by in-medium two-pion and 
kaon exchange processes and denoted by G^^(p). Their contributions to the 
A hyperon self-energy in the medium have been calculated explicitly [32] as 
functions of the nuclear Fermi momentum. This calculation was performed 
using in-medium chiral SU(3) perturbation theory to two-loop order. Following 



^ We retain the notation adopted in Refs. |28|29|30] : space vectors are denoted with 
boldface characters (a;), vectors in isospin space with an arrow (r); Greek symbols 
are used for space-time indices. 

^ The gradient term introduced here is a step beyond our previous hypernuclear 
study [30] which employed a simpler LDA model. 
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the procedures described in Ref. [29] one finds: 

G^^M = AG^ + g,p^+ g, , ^G^ = g^-g^k. (4) 



The hyperon self-energy is derived and calculated in the non-relativistic limit 
at which Lorentz scalar and vector contributions are indistinguishable. Our 
convention is that p^[p) acts with equal share in both scalar and vector 
channels so that its total contribution to the A self-energy is 2G^^(p)p (ig- 
noring the small difference between scalar and baryon densities at this point). 
The density-independent piece AG^ (with constants gi and (72) is associated 
with the regularization of divergent parts of two-pion exchange loop integrals, 
where A ~ 0.7 GeV is a typical cutoff scale. This high-momentum (or short- 
distance) piece is equivalent to a AA^ contact term encoding unresolved short- 
range dynamics at momentum scales large compared to the nuclear Fermi 
momentum. The density dependent terms proportional to (73 and g^ reflect 
the action of the Pauli principle on intermediate nucleons participating in the 
two-pion and kaon exchange processes. The constants gi are deduced from 
[32] , with the following values: 

gi = 2.51 fm^, ^2 = 0.83 fm^, g^ = -0.44 fm^, g^ = 0.84 fml 

Notably, in the terms with non-trivial density dependence representing in- 
medium chiral two-pion and kaon exchange dynamics, g^ and g^ enter with 
alternating signs. 

A second distinct set of contact terms, with coupling constants denoted as 
Gy^^^ and Gg^^\ is introduced to account for the strong Lorentz scalar and 
vector fields of about equal magnitude but with opposite signs, the ones that 
figure prominently in relativistic mean-field phenomenology [M]. In the con- 
text of in- medium QCD sum rules [32], these terms can be associated with 
the leading density dependence of the chiral (quark) condensate, {qq), and the 
quark density (g^g). 

Altogether, the scalar and vector A-nuclear interaction strengths are given as: 
GUp) = gT + Gt,Kip) , G^p) = + Gt^M ■ (5) 



Note that with a typical cut-off parameter A ~ 0.7 GeV, the contact term 
from chiral 27i exchange, AG^ ~ —0.5 fm^, is much smaller than what is 
expected for the individual magnitudes of the scalar and vector contact cou- 
plings Gg^y. For nucleons, Gy ~ ~G^g^ ~ 10 fm^. For the A hyperon we 

also anticipate Gy^^^ ~ —Gg^^\ but with a typical strength of only about half 
that for nucleons. The approximate cancellation in the sum S5 -|- of the 
scalar and vector mean fields produced by G^^v contrasted by their coherent 
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enhancement in the difference Eg — that contributes prominently to the 
spin-orbit couphng. 

In summary, we have the low-density expansion 

Gs,vip) = Giv + 9sP~' +94p'^ , (6) 

with two contact terms, or low-energy constants 

Ssy = G^sf + , (7) 

representing short- distance and "vacuum" dynamics. Apart from the strength 
of the additional surface derivative term, these two constants are effectively 
the only adjustable parameters of the model. Their arrangement in the form 
([7]) is such that individually large "vacuum" parts G^^y are separated from 
the smaller piece AG^ arising from in-medium chiral perturbation theory in 
the non-relativistic limit. Such a separation is not necessary in principle but 
useful in practice. In particular, it helps interpreting the physics content of 
Gg^y in relation to in-medium QCD sum rules, as will be discussed in Section 

[231 

Concerning the coupling strength Dg of the surface (gradient) term in Eq. 
([3]) we follow a similar procedure as in Ref. [29]: Dg will first be treated as 
an adjustable parameter, and its resulting value will then be compared with 
an estimate based on in-medium chiral perturbation theory (cf. Appendix R|) . 
The coupling Dg of the gradient terno could be density dependent as well, 
but it will turn out in practical applications that it is not necessary to go 
beyond the simplest approximation with a constant coupling. 



2.3 Single-particle equations 

The minimization of the hypernuclear ground-state energy leads to a set of 
coupled relativistic equations of Kohn-Sham type for the nucleons and the 
single A-hyperon (see Ref. [30] for further details): 

i^iir) = eiiji{r) withi = n,p,A, (8) 



In non-relativistic calculations such gradient terms have usually been omitted 
on the grounds that finite-range effects are included in an approximate way by the 
use of the empirical charge densities |15j . or they are thought to be absorbed in 
phenomenological non- linear terms in powers of p |15|16|17] . 



-Z7 ■ V + Mi + 70 SV + 
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where ipa{r) are the Dirac wave functions of single-particle orbitals a with 
energies ea for the nucleons and the A hyperon, Mj are their correspond- 
ing masses, and Sy^ denote the vector and scalar self-energies, respectively. 
The hypernuclear ground state is determined by the self-consistent solution of 
the single-particle Kohn-Sham equations ([8]) for a given number of protons, 
neutrons, and for the single A-hyperon. Self-consistency here means that the 
self-energies are functionals of the ground state density calculated in the no- 
sea approximation [Mf36] from the single-particle solutions of the Kohn-Sham 
equations for the nucleons. Because of the explicit density dependence of the 
couplings, rearrangement contributions appear in the vector self-energies 
of protons and neutrons. The A self-energies 

e^ = gUp)p, ^I={gUp) + dIv')ps, (9) 

are then expressed in terms of the vector and scalar ground-state local den- 
sities, p(r) and psi^), of the nuclear core. These densities are determined 
self-consistently together with the wave functions ipa'^i''^) ^ind ip^{r). 

Note that the upper components of the A's Dirac wave function in a given 
orbital, ip^, experience a self-consistent potential = + S^, the sum of 

the vector and scalar self-energies. Given that Gy^'^^ and almost cancel, 

we see that ~ 2G^^(p) p {+ surface term) is of genuine non-relativistic 
origin. Moreover, the dominant part of this average potential comes from the 
short distance (Hartree) piece proportional to AG^, while the sum of the 
and terms gives only a small correction. The lower components of ip^, on 
the other hand, involve the large difference of vector and scalar self-energies, 
Tiy — E;^ ~ {Gy^'^^ — G^^^^) p, that enters in the discussion of the A- nuclear 
spin-orbit coupling, see subsection 14.21 



2.4 Guidance from in-medium QCD sum rules 

The QCD ground state (vacuum) is characterized by condensates of quark- 
antiquark pairs and gluons, an entirely non-perturbative phenomenon. The 
quark condensate (qq), i.e. the vacuum expectation value of the scalar quark 
density, plays a particularly important role as an order parameter of sponta- 
neously broken chiral symmetry. At a renormalization scale of about 1 GeV the 
chiral vacuum condensate is {qq)o — —(240 MeV)'^ ~ —1.8 fm~^ [3H]- Hadrons, 
as well as nuclei, are excitations built on this condensed QCD ground state. 
The density- dependent changes of the condensate structure in the presence 
of baryonic matter are a source of strong scalar and vector mean fields ex- 
perienced by nucleons (and hyperons). In- medium QCD sum rules [35] relate 
the leading changes of the scalar quark condensate {qq)p and quark density 
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(q^q) = 3p, at finite baryon density p, to the scalar and vector self-energies of 
a nucleon (or hyperon) in the nuclear medium. 



The strength of the chiral condensate at normal nuclear matter density, po — 

0. 16 fm~^, is reduced by about one third from its vacuum value. The detailed 
density dependence of this condensate has recently been studied [39] using in- 
medium chiral perturbation theory to three-loop order in the energy density, 

1. e. at a level consistent with the approach employed in the present work. It 
is found that at densities p < po, the leading linear p dependence of {qq)p 
dominates whereas non-linear effects become increasingly important at higher 
densities. Assuming that the nucleon mass in vacuum scales roughly with 
the vacuum chiral condensate, the in-medium scalar and vector self-energies 
of the nucleon can be expressed as 



mi f i mi ri 



Here = {N\mg qq\N) is the nucleon sigma term (~ 50 MeV), is the pion 
mass (138 MeV), and /^r = 92.4 MeV is the pion decay constant. For the quark 
masses we take rriu + fnfi ~ 12 MeV (again at a renormalization scale of about 
1 GeV). The resulting and Sy^°^ are individually large, 300 - 400 MeV 

in magnitude. Their ratio '•°^/Ey'-°'' ~ — (jAr/4(mu + mt^) ~ —1 suggests the 
already mentioned strong cancellation of scalar and vector potentials in the 
single-nucleon Dirac equation. 



The constraints implied by Eq. ffTOj) are admittedly not very accurate, given 
corrections from condensates of higher dimension and uncertainties in the 
values of ctjv and 112^ + 112^. The estimated error for the ratio E^^/S^''' ~ —1 is 
about 20%. Nonetheless, Eq. flTU]) is useful for first orientation when estimating 
the contact couplings G^gy: 

gP ^ ^ -7^^^ ~ -11 fm^ (11) 

m2/2 4(m, + md) 



using (Tat — 50 MeV and + ma ^ 12 MeV. This estimate is actually in 
remarkable agreement with the values determined from a best fit analysis of 
ground-state properties of finite nuclei throughout the nuclear chart {G^ = 
— 11.5 fm^ and Gy^ = 11.0 fm^) [21], the values we use as basic input in the 
present work as well. 

In the case of a hyperon in the nuclear medium, finite-density QCD sum rules 
[ID] predict reduced scalar and vector self-energies of the A while maintaining 
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j]Mo) ^ — Ey*-"-*, though with large uncertainties. We introduce the ansatz 
G'sP = CGfy, (12) 



with a parameter C < 1 controlhng the reduction of the A couphngs relative 
to those for the nucleon. Assuming that only non-strange quarks contribute 
to the contact interactions generated by the condensate background, a naive 
quark model estimate gives ( = 2/3. We leave room for an optimization \3of^ 
of the parameter ( by comparison with empirical A single-particle energies in 
hypernuclei. In fact, while the quark model value is an option, the result of 
the detailed QCD sum rule analysis [ID] is closer to C = 0.4. 



3 Single particle states of the A 



Given the explicit input for the long and intermediate range kaon and two- 
pion exchange interactions, the only remaining unknowns are the strengths of 
the short-distance (contact) terms, G^^y, and of the surface gradient term, 
Dg. These will now be fixed by detailed fits to the empirical single particle 
orbits of the A in selected hypernuclei. The contact terms are specified by the 
parameter ( relating the scalar and vector mean fields of the A to those of the 
nucleons, see Eq. flT^ . The remaining AG^ involves the cutoff scale A. One 
expects this cutoff scale to be around 0.7 GeV, subject to possible further 
fine-tuning. 

We proceed as follows. The values of A and are adjusted for different values 
of ( (using 0.4,0.5 and 2/3, i.e. ranging from the QCD sum-rule estimate to 
the naive quark- model prediction), by performing a least-squares fit to the 
empirical single-A energy levels: 

( th _ex'p \ 2 



Here e„ and e'^^ are the theoretical and experimental single-A energies, re- 
spectively, with uncertainties 6e^^. For the set of experimental energies in 
Eq. ([13]) we choose the s- and p- levels of ^^O and ^ofPb (see Table |2]|3. An 
uncertainty of 3% is assumed for the s-states, and 5% for the j)-states. For 

^ In Ref. [30] we have used x instead of C for the ratios G^^^^ /G^'^\ 
^ In a recent Jlab Hall-A report [H], Cusanno et al. studied the single particle 
spectrum of ^^N and found for the Is binding energy the value 13.76 it 0.16 MeV, 
in apparent disagreement with the corresponding value for ^^O if charge symmetry 
would be approximately realized. 
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Table 1 

Best-fit values of the contact term AG^ (determined by the cut-off scale A) and 
of the strength of the surface (gradient) term, for different choices of the ratio 
( = G^'^y^/G^^^y (input: G^ = -11.5 fm^, gJ?^ = 11.0 fm^ from Ref.[29]). The fits 
are performed by reproducing single particle energies of s- and p-states in ^^O and 
^'^^Pb. Also shown is the resulting depth of the A central potential. 



c 


AG^ (fm2) 


A (MeV) 


Dl (fm^) 


(MeV) 




0.4 


-0.55 


718.6 


-0.304 


-34.7 


3.41 


0.5 


-0.56 


721.0 


-0.340 


-35.4 


4.34 


2 
3 


-0.58 


726.1 


-0.415 


-36.6 


2.86 



the nucleoli sector the FKVW energy density functional [22] is used, without 
readjustments of any parameters. The fits are summarized in Table [H where 
we display the optimal values of the cut-off scale A (or equivalently, the con- 
tact term proportional to AG^), and the strength of the gradient term Df , for 
C = 0.4, 0.5 and 2/3. The corresponding values are given in the column on 
the right. Note that at this stage the data set used in the fit does not include 
the (small) splittings between the p-shell spin-orbit partner states, but only 
the average energies of the p-orbitals. 

Quite acceptable fits are obtained with AG^ in the range —(0.55 — 0.58) 
fm^, producing a central A single particle potential U^entrai = 2AG^p(r) ^ 
— (35 — 37) MeV (p(r)/po)- The required cutoff scales A are close to A ~ 0.7 
GeV [32], almost independent of the ratio C, = Gg^y /G^g\/. In view of the 

cancellation between scalar and vector contact terms G^g^ and Gy\ this is of 
course not surprising. 

The inclusion of the surface term proportional to VpA ■ Vp turns out to be 
important. An attempt to fit the A energy levels without this term results in 
an unacceptably large value (^ 200), because it is not possible to simulta- 
neously reproduce levels in light and heavy hypernuclei. For the strength of 
this term, Dg, one could directly use the theoretical estimate —0.56 fm^ from 
in-medium chiral perturbation theory (cf. Appendix . The resulting agree- 
ment with data is then significantly improved, but not yet optimal [x^ ^ 10). 
In particular, one still encounters overbinding for heavy systems. The best fit 
values, ~ -(0.3 - 0.4) fm^ are nevertheless still in remarkable qualitative 
agreement with chiral theory. 

In this work all calculations are performed assuming spherical symmetry and a 
simplified configuration: a closed core of nucleon pairs plus a single A-hyperon. 
These are commonly used approximations and recent investigations have con- 
firmed their validity. The study of Ref. [T7] has shown that deformation effects 
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need to be taken into account only in very light systems which are not the sub- 
ject of the present investigation. Of the hypernuclei considered in this work, 
only ^^C is sensitive to deformation, but recent RMF studies [12] have pointed 
out that the stabilizing effect of the A tends to restore spherical symmetry in 
the hypernuclear system. Effects related to the odd number of nucleons in hy- 
pernuclei appear to be rather small [17] . In the case of hypernuclei with open 
nucleon shells, we include pairing correlations described in the BCS approx- 
imation with empirical pairing gaps [13] (except for ^^C, which is calculated 
in the closed-shell approximation). 



4 A-nuclear spin-orbit coupling 



Having specified the input for the A-hypernuclear central and surface po- 
tentials, we can now concentrate on a more detailed investigation of the A- 
nuclear spin-orbit interaction, one of the central themes of hypernuclear spec- 
troscopy. As already indicated in the introduction, A-hypernuclei feature ex- 
tremely small energy spacings between spin-orbit partner states, as compared 
to the large spin-orbit splittings in ordinary nuclei [T][2|3] . In this section we 
first briefly summarize the available data and review the current theoreti- 
cal approaches. The final subsection presents and discusses results based on 
our novel interpretation of the smallness of the A-nucleus spin-orbit coupling 



4-1 Brief summary of experimental results 



Detailed informations on A spin-orbit splittings derive from measurements of 
light hypernuclei with nucleons in the j>-shell. The BNL 929 experiment |44j . 
using ^^C as target, reported an energy splitting between the pi/2 and P3/2 
A levels Ae^(p) = 152 ± 54 ± 36 keV. At present this experimental result 
provides the most convincing evidence for the small spin-orbit coupling in A- 
hypernuclei. It has been corroborated by experiments on heavier hypernuclei. 
For ^^O a combined analysis of {K~ , 7r~7) and (vr^, K~^) experimental spectra 
(experiments KEK E336 [15] and CERN-SPSII [16]), allowed the determina- 
tion of the very small energy level splitting between the 2"*" and 0^ states: 
Ae(2+ - 0+) = 40 ±320 keV. Motoba et al. [17] used the linear dependence of 
the relation between Ae(2+ — O"*") and the p-\eve\ splitting Ae^{p), to estimate 
Ae^{p) = 300 — 600 keV0. An updated analysis [B] of data from the experi- 
ment ^^Y(7r+, K+flY [18] has confirmed that the spin-orbit interaction is also 

^ A more recent study by Hashimoto et ai, lowers this estimate to: —800 < 
Ae^(p) < 200 keV (cf. Fig. 23 in Ref. [I]). 
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strongly suppressed in heavy hypernuclei: Ae^(/) ^ 200 keV, Ae^{d) ^150 
keV and Ae^{p) ~ 90 keV. So far these are the only empirical spin-orbit 
splittings outside the region of j>-shell nuclei. 



4-2 Previous theoretical studies 



The framework of energy density functionals provides an accurate description 
of hypernuclei over the whole mass table [15] . At the same time light systems 
(e.g. j>-shell hypernuclei) are also successfully decribed using phenomenological 
hyperon-nucleus potentials such as the one introduced by Dalitz and collab- 
orators jlHj. More recent ab-initio calculations start from realistic free YN 
potentials, as for example the Nijmegen studies based on meson-exchange in- 
teractions [SnEI], or the Jiilich potential based on SU{3) chiral perturbation 
theory [52]. However, all these realistic YN potentials tend to strongly overes- 
timate the AA^ spin-orbit interaction in hypernuclei. For instance, recent cal- 
culations based on the cluster model approach [53] predict Ae^(p) ~ 390 — 960 
keV for i3c[7]. 

For first orientation, recall the Walecka-type relativistic mean-field approach 
to nuclei [51] in which short-range dynamics is parametrized by the exchange of 
a phenomenological scalar boson (a) and a vector boson (u) between nucleons. 
The corresponding mean fields produce large scalar (S) and vector {V) nucleon 
self-energies. The effective spin-orbit potential [57] is obtained in the non- 
relativistic limit of the single-particle Dirac equation: 



Vso = ^J-^VUr)]l-s, (14) 




where the large spin-orbit potential Vis arises from the difference of the vector 
and scalar potentials: S^/(r) = g^jUJo^r) and ^^(r) = g„a{r), with Sy(0) ~ 
330 MeV and ^^(0) ~ -400 MeV [55115^11^ . Exphcitly, 

M 

1^, = -— (Sv^-S<j), (15) 



With a recent improvement of the Nijmegen potential, called ESC06 [5l], it ap- 
pears possible to obtain values of Ae^ that are closer to data. However, Nijmegen 
potentials such as ESC04, fail to reproduce the single-A binding energies [55]. Thus 
the microscopic interpretation of the A- nuclear spin-orbit interaction, starting from 
one-boson exchange YN interactions, remains an open problem. For the latest ver- 
sion (ESC07) including also contributions from quark degrees of freedom, see [56] , 
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where Me// is an effective mass specified as [59] 
M,// = M-i(Sv.-S5). 



(16) 



At this point it is useful to clarify the close correspondence between a Walecka 
type phenomenology (or an in-medium QCD sum rule approach) of strong 
scalar-vector mean fields, and certain contact terms appearing at next-to- 
leading order (NLO) in the effective field theory description of the NN inter- 
action. A combination of NLO contact terms (with two derivatives, see Eq. 
(5) in Ref. [SI]) generates a Galilei invariant spin-orbit interaction 

— Cs (cTi + 0-2) ■ 



P2) X {P1-P2) 



(17) 



where C5 is given by a linear combination of P-wave low-energy constants: 

C5 = -^[2C(3Po) + 3C(3Pi) - 5C(=^P2)] . (1 
167r 



It has been shown in Ref. [61] that the strength of the spin-orbit term derived 
from the free NN interaction agrees quantitatively with the one required in 
nuclear shell-model calculations. It has furthermore been demonstrated in Ref. 
[62] that the spin-orbit contact term (ITS!) . in a relativistic Dirac-Briickner 
(DBHF) calculation, generates the strong scalar and vector mean fields in 
the combination S5 — Sy. Thus the inclusion of contact terms representing 
condensate background fields does not at all spoil the consistency of in-medium 
chiral calculations based on the LO vrA^ Lagrangian. These background fields 
are just equivalent reflections of short distance A^A^ dynamics. 

The scalar-vector approach can be transcribed analogously for hypernuclei. 
Empirical spin-orbit splittings are reproduced by simply assuming much weaker 
couplings between the A and the exchanged bosons. In particular, in one of 
the early studies of this type [7], a reduction of 1/3 was suggested for the A 
potentials with respect to the corresponding nucleon self-energies 

^s = l^s and S^ = ^Sy. (19) 



In contrast, the naive quark model assumes that the non-strange quarks cou- 
ple to the a and u mean-fields (the s-quark spectator hypothesis) and suggests 
a reduction factor of 2/3. With this value, however, it is not possible to repro- 
duce the empirical spin-orbit splittings in A-hypernuclei. A possible solution 
proposed in Ref. [8] involved an additional strong tensor coupling term in the 
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uoK interaction Lagrangian 



2Ma 



i)A(y^'''ipAdyOJ^ . 



(20) 



This additional term modifies the effective A spin-orbit potential as follows: 



For f^/gt = -1 the potential V^,(A) = (2/^/^?;^ + 1)S^ - is now very 
small compared to that for the nucleon (see Ref. [TU] for more details, and 
Figs. 2 and 3 of Ref. [H]). 

While phenomenological studies, based on the assumption of a strong uA. ten- 
sor coupling [11], worked in reproducing the empirical single- A levels for a 
number of hypernuclei, they did not offer a consistent microscopic explana- 
tion for the spin-orbit suppression in A-hypernuclei. Realistic YN potentials 
(NSC97) that reproduce phase shift data suggest significantly weaker uK ten- 
sor couplings [50]. The older Nijmegen D and F potentials, for example, give 
-0.12 and -0.54 for the ratio fl^/g^ [9], respectively, considerably smaller in 
magnitude than the value — 1 required in Eq. (l2T|) to reproduce the empirical 
spin-orbit splittings. 

Alternative microscopic models of hypernuclear spectroscopy have also been 
developed more recently. A synthesis of quark-model and relativistic one-boson 
exchange picture has been established by the Quark-Meson Coupling (QMC) 
model [T2p3] . In this model Ko,a arises entirely from Thomas precession. To 
obtain the correct spin-orbit splittings, a piece —j^^g^-^uj{r)l ■ s must be 
included in addition to the self-consistent calculation of single-A energy levels. 
Lenske et al. [H] have developed a density-dependent relativistic framework in 
which the A-meson couplings are partly determined from a AN T-matrix, and 
partly fitted to a selected set of data. In that approach the spin-orbit energy 
splittings Ae^ display a uniform dependence on the nuclear mass number A 
(cf. Fig. 7 of Ref. [S]). This is in contrast to the QMC results, where Ae^ ~ 
for all hypernuclei and any A orbital [12]. The spin-orbit energy splittings 
computed in [H] still overestimate the empirical values, and this led those 
authors to raise the quest for an additional reduction mechanism of the A- 
nuclear spin-orbit force. 




(21) 
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4-3 The A-nuclear spin-orbit interaction from chiral SU(3) two-pion exchange 



Starting from in-medium chiral SU (3) dynamics for the A hyperon in nuclear 
matter and its translation to a hypernuclear energy density functional 

[3U] . we have identified two basic competing mechanisms at the origin of the 
unusually small A-nuclear spin-orbit interaction: strong scalar-vector mean 
fields acting coherently, and the spin-orbit force of opposite sign [32] induced 
by the second order pion exchange tensor interaction with an intermediate 
S hyperon. A third prominent contributor to the spin-orbit force in ordinary 
nuclei, namely the three-body interaction of Fujita-Miyazawa type, has no 
counterpart in single-A hypernuclei, as already mentioned. 

Let us briefly recall the steps leading to the "wrong sign" spin-orbit interaction 
from two-pion exchange [32]. Consider a A that scatters in slightly inhomoge- 
neous nuclear matter from an initial momentum p — q/2 to a. final momentum 
p + q/2. One identifies a spin-orbit term T,f^{kp) = ^U^i^kp) cr ■ {q x p), in 
the spin-dependent self-energy of the A, with 



' D9A 



\pi\<kF 



d^Pi 
(2^ 



\P2\>kF 
iPl 



d^P2 



P2] 



B 



[ml + (pi - P2YY [^^ + P2 ^ - Pi ■ P2] 



.(22) 



It depends only on known SU (3) axial-vector coupling constants [D = 0.84, 
= 1.3) and on = (Ms — Ma)Mb which involves the (small) mass difference 
between the S and the A hyperon. The average baryon mass Mb = 1.05 GeV 
appearing in the numerator is a reminder that this spin-orbit interaction has 
a non-relativistic origin. The momentum space loop integral (122|) is finite and 
hence model independent, in the sense that no regularizing cutoff is required. 
The spin-orbit coupling strength at saturation density, U^{kp^) ^ —15.1 MeV 
fm^ at kp^ ~ 1.36 fm~^, has a sign opposite to the standard spin-orbit cou- 
pling strength. Evidently, this contribution to the A spin-orbit potential tends 
to cancel the contribution from the strong scalar- vector mean fields. However, 
contrary to the case of nucleons in ordinary nuclei (see Fig. [1]), this "wrong 
sign" spin-orbit interaction is not compensated in turn by a three-body spin- 
orbit interaction [52|55] . Thus the smallness of the A- nucleus spin-orbit finds 
its natural explanation in terms of an almost complete cancellation between 
short-range background mean-field contributions and longer range terms gen- 
erated by 27r-exchange. 

One might argue that the ?7/^ of Eq. (!22|) . derived for infinite nuclear matter 
with its continuous single particle energy spectrum, overestimates the mag- 
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Fig. 1. Schematic picture of contributions to the spin-orbit interaction of a nucleon 
in nuclei and a A in hypernuclei. First row: contact terms summarizing strong scalar 
and vector fields and unresolved short-distance dynamics. Second row: iterated pion 
exchange (mostly second order tensor interactions) with Pauli-blocked intermediate 
nucleon states. Third row: three body terms generated by 27r exchange with interme- 
diate excitation of a virtual A isobar and in-medium insertions (resembling Fujita 
and Mijazawa |31j). This mechanism has no counterpart for A hyperons since there 
is no filled hyperon Fermi sea. Signs of the different contributions are indicated. 



nitude of the two-pion exchange spin-orbit force considerably. This turns out 
not to be the case. In order to examine this question, consider instead a finite 
cubic box of length 2L with a discrete momentum spectrum, p = / L)n and 
n e Z^. Momentum space integrals are now replaced by sums, / d^p/ {2n)^ 
(1/8L3) J2n- The density p is 1/{2L^) times the sum over all occupied, discrete 
single particle states. Now take, for example, a finite system with A = 48 for 
which L ~ 3.35 fm. The result is UH ~ —14.7 MeV fm^, to be compared with 
-15.1 MeV fm^ for nuclear matter. 

The quantitative effect of in-medium two-pion exchange on the energy spac- 
ings between single-A spin-orbit partner states, is calculated in first-order 
perturbation theory: 



Ae^ = {a\Ant\a) , 



(23) 
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for each hyperon orbit a, with 



^nt = ^^^^^i-\'-^\cT.l. (24) 




UtAkf) d (p{r 



p(o) 



The A single particle states |a) are determined by self-consistent solutions of 
the Dirac equation ([8]), and p(r) is the corresponding ground-state nuclear 
density, also computed self- consistently, with p^^^ = p(r = 0). 
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Fig. 2. Calculated A single particle energy levels in ^O. Left panel shows the spin-or- 
bit splitting from short-distance dynamics (scalar and vector mean fields) only; right 
panel: after inclusion of the spin-orbit term from in-medium two-pion exchange. 

A demonstration of the counterbalance between short-distance contact term 
and intermediate range two-pion exchange contributions to the spin-orbit split- 
ting is shown in Fig. [2] for the example of ^j^O. It shows how the spin-orbit term 
from the in-medium two-pion exchange (second-order tensor) force basically 
compensates the contribution from scalar-vector mean fields. 

We close this subsection with a note on corrections to the chiral spin-orbit 
potential from the 5't/(3) decuplet. A study along these lines has been per- 
formed in Ref. [65] which extended the work of Ref. [32] and considered also 
A and S* intermediate states in the two-pion exchange processes. Corrections 
from virtual excitations of such states, despite their large masses, turned out 
to be non-negligible because of their strong couplings to the pions and the 
baryons octet. The results found previously in Ref. [32] were basically con- 
firmed: the spin-orbit potential from two-pion exchange is a non-relativistic 
effect as demonstrated by its proportionality to the baryon mass. It is model 
independent as it does not require regularization and depends only on physical 
parameters (masses and coupling constants). We have examined this extended 
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[//g in our approach as well and briefly report results at the end of the next 
section. 



5 Results 



The previously decribed framework, with parameters fixed once and for all by 
the s- and p-states of ^^O and Pb, is now tested systematically in comparison 
with the empirical energies for all known single-A levels of six hypernuclei, 
from light to heavy. The input from Tabled] with ( = 0.5 is used ^ I. together 
with the FKVW parameters [22] in the nucleon sector of the energy density 
functional but, as will shortly be shown, the results are not very sensitive to 
a particular choice of one of the three parameter sets for the AA^ couplings. 

Our results are summarized in Tables [2] and [3] and compared with those of six 
different calculations: Quark Meson Coupling (QMC) [T2IT3] . Fermi Hyper- 
netted Chain (FHNC) [18], Skyrme (SK) [16], Brueckner-Hartree-Fock (BHF) 
[TU] with the Njimegen SC97F potential [SD], relativistic mean field models 
with a tensor coupling [TT] (RMFI with f^/g^ = —1) and density-dependent 
interactions [11] (RMFII). The single particle energy spectra are calculated in 
the approximation of a closed even-even nucleon core -|- A-hyperon, i.e. the 
theoretical spectra0 correspond to ^^C, ^O, ^Ca, AZr, ^^^Ce and ^°^Pb. 

Comparing our (FKVW) results to those of a relativistic mean-field calcula- 
tion (RMFI) that uses a strong A-nuclear tensor force with fl^ /g^ = —1, one 
observes differences in the spin-orbit splittings of A orbitals with large angu- 
lar momentum /. For example, the FKVW spin-orbit splitting, though very 
small on an absolute scale and well within experimental errors, tends to be 
more than twice the RMFI splitting for a A in /- and gf-orbitals. The mech- 
anisms at work in FKVW and RMFI involve of course very different physics. 
A technical difference arises because the FKVW spin-orbit splitting is treated 
perturbatively, whereas the RMF calculations are performed self-consistently. 



We recall that the results are not sensitive to a particular choice of parameter 
sets in Table [TJ 

^ for the calculation should not be considered fully realistic due to the closed 
shell approximation. 
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Table 2 

Binding energies (in MeV) of single-A levels in ^^C, ""^^O, ^^Ca and *^Y. Experimen- 
tal energies [1] are shown in comparison with the results of the present calculations, 
using the input parameters of Table[T]and C = 0.5 (column FKVW). Also listed are 
results of five different models: Quark Meson Coupling (QMC) |12|13| . Fermi Hyper- 
netted Chain (FHNC) [E], Skyrme (SK) [H], Brueckner-Hartree-Fock (BHF) [H] 
with the Njimegen SC97F potential [50], and RMF models with a tensor coupling 
[TT] (RMFI with f^/g!^ = -1) and density-dependent couphngs [U] (RMFII). 



AT 1 

iNucieus 




Expt. 


FKVW 


QMC 


FHNC 


SK 


BHF 


RMFI 


T~» TV /TT^TT 

RMFII 


13 p 


lSl/2 


11.38 ±0.05 


12.3 


- 


8.3 


11.7 


13.7 


12.5 


11.7 




1^3/2 


0.38 ±0.1 


0.1 






0.9 


1.4 


1.1 


1.1 




lPl/2 




0.0 










0.8 


0.0 


16 Q 

A^ 


lSl/2 


12.42 ±0.05 


12.6 


16.2 


12.00 


13.3 


15.5 


12.9 


12.8 




1^3/2 


1.85 ± 0.06 


2.0 


6.4 


1.8 


3.0 


3.7 


3.3 


2.8 




lPl/2 




1.9 


6.4 








3.0 


1.4 


^°Ca 


lSl/2 


20.0 ± 1.0 


18.9 


20.6 


20.0 


18.0 


20.7 


19.0 


17.6 




1^3/2 


12.0 ± 1.0 


10.1 


13.9 


10.6 


10.1 


11.5 


10.7 


9.1 








10.1 


13.9 








10.5 


7.8 




1^5/2 


1.0 ± 1.0 


1.6 


5.5 


1.6 


1.6 


2.0 


2.7 


1.5 




lf^3/2 




0.9 


5.5 








2.4 


1.5 


89 V 
A ^ 


lSl/2 


23.1 ±0.5 


23.4 


24.0 


23.3 


21.1 


24.1 


23.7 


23.2 




lP3/2 


16.5 ±4.1 


17.2 


19.4 


16.9 


15.6 


17.8 


17.6 


17.2 




lPl/2 




17.2 


19.4 








17.4 


16.3 




1^5/2 


9.1 ± 1.3 


10.2 


13.4 


10.1 


9.1 


10.4 


10.7 


10.3 




1^^3/2 




9.8 


13.4 








10.5 


8.9 




1/7/2 


2.3 ± 1.2 


2.8 


6.5 




2.1 


2.4 


3.7 


3.1 




1/5/2 




2.0 


6.4 








3.4 


1.0 



The present calculations (column FKVW) are evidently in very good agree- 
ment with data, and comparable in quality or superior (especially for heav- 
ier hypernuclei) to the energy spectra calculated with other approaches. The 
QMC model in its original form included a phenomenogical spin-orbit correc- 
tion and the Pauli-blocking effect at the quark level. Without these corrections 
the resulting energy levels show a strong overbinding (cf. Tab. 4 in Ref. |12j). 
A very recent improvement [13] solved the overbinding problem, introducing 
the scalar polarizability of the nucleon in a self-consistent way instead of the 
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Table 3 

Binding energies (in MeV) of single- A levels in ^■^^La and ^^f Pb (continued from 
Table E]). 



Nucleus 




Expt. 


FKVW 


QMC 


FHNC 


SK 


BHF 


RMFI 


RMFII 






24.5 it 1.2 


24.7 






22.1 


25.3 


25.2 


25.2 




lP3/2 


20.4 ±0.6 


20.0 


— 


— 


17.9 


20.5 


20.4 


20.5 




lPl/2 




20.0 










20.4 


20.2 




lf^5/2 


14.3 ±0.6 


14.3 








12.8 


14.5 


14.8 


14.9 




1^3/2 




14.1 










14.6 


14.1 




1 f 

J^/7/2 


8.0 ±0.6 


Q n 
o.U 


— 


— 


6.9 


7.8 


o.D 


0.0 




1/5/2 




7.4 










8.4 


7.1 




109/2 


1.5 ±0.6 


L.O 


— 


— 


0.6 


0.6 


Z.4 


z.z 




I57/2 




0.5 










2.0 


0.2 


Pb 


lSl/2 


26.3 ± 0.8 


25.8 


26.9 


27.6 


23.1 


26.5 


26.5 


27.2 




lP3/2 


21.9 ±0.6 


22.0 


24.0 


22.8 


19.6 


22.4 


22.7 


23.4 




lPl/2 




22.0 


24.0 








22.6 


23.1 




1^5/2 


16.8 ±0.7 


17.4 


20.1 


17.4 


15.4 


17.5 


18.0 


18.5 




1^3/2 




17.3 


20.1 








17.9 


17.9 




1/7/2 


11.7 ±0.6 


12.2 


15.4 




10.5 


11.8 


12.7 


13.2 




1/5/2 




11.8 


15.4 








12.5 


12.1 




I99/2 


6.6 ±0.6 


6.5 


10.1 




5.1 


5.6 


7.1 


7.5 




I57/2 




5.8 


10.1 








6.9 


5.8 



Pauli blocking correction. In Tables [2] and [3] we have included the latest update 
of these calculations. 



Figure [3] provides a further test of the sensitivity of calculated single-A ener- 
gies with respect to a variation of the ratio ( = G^^y /G^g\r between contact 
terms representing the in-medium condensate background fields for the hy- 
peron and the nucleons. For the six hypernuclei listed in Tables |2] and [3], the A 
binding energies calculated with the FKVW parameters plus the three best-fit 
parameter sets from Table [T] that determine the AA^ couplings, are plotted as 
functions of the mass number and compared with empirical energies. Calcula- 
tions with all three parameter sets of the AA^ interaction reproduce the data 
with high accuracy for a wide range of hypernuclear masses, using a reason- 
able band width of ( values between the QCD sum-rule estimate {( ~ 0.4) 
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Fig. 3. Binding energies of the A in different (s,p, ...) orbitals of six hypernuclei 
(cf. Tables [2] and [3]) , calculated with the FKVW density functional using the three 
parameter sets for the AN couplings (cf. Table [I]) . Results are plotted as functions 
of the mass number and compared with experimental energies [1]. Also shown is a 
Woods-Saxon fit [15] (dashed curves) to guide the eye. 



and the naive quark- model = 2/3). 



In Tab. HI we show the theoretical spin-orbit splittings for A p-levels in compar- 
ison with empirical values (second column) and other relativistic calculations 
with (fourth column) and without (fifth column) tensor coupling. 



Concerning the A-nuclear spin-orbit splittings, these calculations have all been 
performed using the chiral SU(3) two-pion exchange AA^ interaction with just 
the E hyperon in the intermediate state. Adding terms with decuplet interme- 
diate states as studied in Ref. [65] do not change the picture in any significant 
way, except that, looking at Ae^(p), there is then a preference for the smaller 
ratio ( ~ 0.4, consistent with the in-medium QCD sum rule analysis. 
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Table 4 

P-shell spin-orbit splittings A = Ae^{p) for six hypernuclei ""^^O, ^Ca, ^^Y, 

^^La, ^5fPb). Experimental values [H], or empirical estimates |1|47I48| . are shown 
in comparison with our theoretical predictions (FKVW), using a broad range of 
parameters (see Eq. [T2|) . and other relativistic calculations with (RMFI [11]) or 
without (RMFII [H]) tensor coupling. All energies are given in keV. The asterisk 
means that a local fit has been necessary. 



Nucleus 


Exp. 

A [keV] 


FKVW 
(0.4 < C < 0.66) 


RMFI [H] 


RMFII m] 


13 p 


152 ± 54 ± 36 [44j 


-160 < A < 510 


310 


~ 1100* 


10 


300 < A < 600 [47] 
-800 < A < 200 pj 


-210 < A < 490 


270 


~ 1400 


ICa 




-140 < A < 420 


210 


~ 1400 


89 V 
A ^ 


90 m] 


-40 < A < 180 


110 


^ 700 


ifLa 




-20 < A < 80 


50 


^ 300 


208pb 




-20 < A < 70 


50 


- 300 



6 Summary and conclusions 

1. A previously derived (FKVW) relativistic nuclear energy density functional, 
with constraints from low-energy QCD, has been generalized to hypernuclei. 
In-medium chiral SU(3) dynamics is implemented at three-loop order in the 
energy density, with explicit treatment of two-pion (and kaon) exchange AN 
interactions in the presence of the filled Fermi sea of nucleons. 

2. Strong scalar and vector fields experienced by the A are considered to 
emerge from its coupling to density-dependent quark condensates constrained 
by in-medium QCD sum rules. They are manifest in the form of contact terms 
that can be related to low-energy constants in a (chiral) effective field theory. 
The corresponding scalar and vector Hartree potentials have opposite signs 
and cancel in the central mean field. Their magnitudes are about 0.4 - 0.5 
of the corresponding mean fields for nucleons in nuclei. These reduced mean 
fields are consistent with the (admittedly uncertain) QCD sum rule analysis 
for nucleons and hyperons in a nuclear medium. 

3. A A-nuclear surface coupling term, that appears in the gradient expansion 
of a density functional for finite systems, is generated model-independently 
from in-medium chiral SU(3) perturbation theory at the two-pion exchange 
level. This term proves to be important in obtaining good overall agreement 
with A single-particle spectra throughout the hypernuclear mass table. 
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4. The chiral two-pion exchange AN interaction in the presence of the nuclear 
core generates a (genuinely non-relativistic, mo del- independent) contribution 
to the A-nuclear spin-orbit force. This longer range contribution counterbal- 
ances the short-distance spin-orbit terms that emerge from scalar and vec- 
tor mean fields, in just such a way that the resulting spin-orbit splitting of 
A single particle orbits is extremely small. A three-body spin-orbit term of 
Fujita-Miyazawa type that figures prominently in the overall large spin-orbit 
splitting observed in ordinary nuclei, is absent for a A in hypernuclei, simply 
because in a single-A hypernucleus there is no Fermi sea of hyperons. 

5. The confrontation of this highly constrained approach with empirical A 
single-particle spectroscopy turns out to be quantitatively successful, at a level 
of accuracy comparable to that of the best existing hypernuclear many-body 
calculations. The resulting A-nuclear single-particle potential has dominant 
Hartree term with a central depth of about —30 MeV, consistent with earlier 
phenomenology. 

6. While a relativistic framework has been used here for practical convenience, 
this is not mandatory. Given the cancellation of relativistic scalar and vector 
terms in the central mean field, their coherent effect in building up part of the 
spin-orbit force can be translated into derivatives of contact terms at next- 
to- leading order in an equivalent, non-relativistic effective field theory. The 
important compensating two-pion exchange mechanism that renders the over- 
all spin-orbit coupling for the A hyperon so abnormally small, is entirely of 
non-relativistic origin, as well as the similarly important surface gradient term 
that adds to the A-nuclear potential. 
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A Appendix: Gradient term of the A-nuclear density functional 

We briefly outline the chiral EFT estimate for the strength parameter Dg of 
the gradient term in the interaction part of the A-nuclear density functional 
Eq. m- 
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q/2 



-q/2 




l + q/2 



Pi + q/2 



Fig. A.l. Two-pion exchange between a A-hyperon and nucleons, with a S-hyperon 
and Pauh blocking effects in the intermediate state. The horizontal double lines 
represent the filled Fermi sea of nucleons. 



This parameter can be evaluated from the spin-independent part of the self- 
energy of a A-hyperon interacting with weakly inhomogeneous nuclear matter. 
Let the A scatter from initial momentum —q/2 to final momentum +q/2, as 
shown in the two-pion exchange diagram of Fig. lA.li For small momentum 
transfer q to the inhomogeneous nuclear medium, the A mean-field potential is 
modified by a correction term proportional to q^. In the energy density func- 
tional the factor q^ generates (via Fourier transform) the product of density 
gradients: VpA ■ Vp (with p = 2k\/?>n'^ the nuclear density). 

Consider now the two-pion exchange between the A and the nucleons in the 
medium. The small momentum transfer q enters both the pion-baryon vertices, 
and the pion propagators. Expanding the diagrammatic expression for the A 
self-energy up to order g^, one finds: 



^) ^^_D''g\MB f d^pidH P {41^ + 91^ ml + 3 mi] 



6/4 J (27r)6 (m2 +/2)4[A2 + /2_^.p^] 

\pi\<kp 



+ Pauli blocking terms , (A.l) 



with the axial vector coupling constants D = 0.84 and qa = 1-3. Here 
Mb = 1047 MeV denotes an average baryon mass, and the SA mass-splitting 
Ms — Ma = A'^/Mb = 77.5 MeV has been rewritten in terms of the small 
scale parameter A = 285 MeV. The function D^^"^ (p) determines the (possibly 
density-dependent) strength of the gradient term. The Pauli blocking correc- 
tion in Eq. (lA.ll) is obtained by reversing the sign, and substituting the loop 
momentum as I = pi — p2 with \p2\ < kp- In this way the Pauli-blocked inter- 
mediate nucleon states in the filled Fermi sea get properly removed from the 
loop integral. Note also that the loop integral in Eq. (lA.ip is finite as it stands 
and does not require any regularization. It could even be solved in terms of 
elementary functions (an arctangent and square-roots). 



The density-matrix expansion of Negele and Vautherin [66] contributes a term 
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proportional to V^p to the in-medium insertion for an inhomogeneous many- 
nucleon system. Taking this term into account one obtains the following ad- 
ditional contribution to the strength function: 



48/^4 J (27r)6 {ml + Py[A^ + P-l-pi] 

\pi\<kF 

+ Pauli blocking terms (A. 2) 

This combined loop and Fermi sphere integral is also convergent as it stands. 
The weighting factor 5pi — 3kp ensures that the divergent constant from the 
loop integral (scaling e.g. with an ultraviolet cutoff) disappears in the final 
result . 

Summing up all four contributions from Eqs. ( lA.ll) and ( ]A.2I) . the resulting 
density-dependent strength function pD^^^ (p) is shown in Fig. IA.2[ Since the 
density dependence is almost linear, the function D^^^ (p) can be approximated 
by its value at nuclear matter saturation density po = 0.16 fm~^ : 

D^I\pq) = (139.3-21.0- 10.6 + 1.3) MeV fm^ = 109 MeV fm^ (A.3) 



Note that this estimate for the strength of the gradient term is a model- 
independent and parameter-free result derived from the long-range 27r-exchange 
with a S-hyperon in the intermediate state. Comparing the gradient term for 
the A-hyperon with the corresponding term in A-nuclear density functional 



26 



Eq. ([3]), we obtain the following relation: 



— ■ 



(A.4) 



From Eq. (lA.SP it then follows: 



= -0.56 fm^ , 



(A.5) 



remarkably close in value to the chiral EFT prediction for the strength of the 
corresponding gradient term in the nucleon sector of the density functional: 
Ds = -0.7 fm^ . 
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